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Counterion adsorption on a flexible polyelectrolyte chain in a spherical cavity is 
considered by taking a "permuted" charge distribution on the chain so that the "ad- 
sorbed" counterions are allowed to move along the backbone. We compute the degree 
of ionization by using self-consistent field theory (SOFT) and compare with the pre- 
viously developed variational theory. Analysis of various contributions to the free 
energy in both theories reveals that the equilibrium degree of ionization is attained 
mainly as an interplay of the adsorption energy of counterions on the backbone, 
the translational entropy of the small ions, and their correlated density fluctuations. 
Degree of ionization computed from SOFT is significantly lower than that from the 
variational formalism. The difference is entirely due to the density fluctuations of 
the small ions in the system, which are accounted for in the variational procedure. 
When these fluctuations are deliberately suppressed in the truncated variational 
procedure, there emerges a remarkable quantitative agreement in the various con- 
tributing factors to the equilibrium degree of ionization, in spite of the fundamental 
differences in the approximations and computational procedures used in these two 
schemes. Furthermore, it is found that the total free energies from the truncated 
variational procedure and the SCFT are in quantitative agreement at low monomer 
densities and differ from each other at higher monomer densities. The disagreement 
at higher monomer densities is due to the inability of the variational calculation 
to accurately compute the solvent entropy at higher concentrations. A comparison 
of electrostatic energies (which are relatively small) reveals that the Debye-Hiickel 
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estimate used in the variational theory is an overestimation of electrostatic energy 
as compared to the Poisson-Boltzmann estimate. Nevertheless, since the significant 
effects from density fluctuations of small ions are not captured by the SCFT, and 
due to the close agreement between SCFT and the other contributing factors in the 
more transparent variational procedure, the latter is a better computational tool for 
obtaining the degree of ionization. 

I. INTRODUCTION 

Counterion adsorption in polyelectrolyte solutions is one of the fundamen- 
tal problems in polyelectrolyte physics, which has been a topic of extensive 
researcW-^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ for decades. Historicallyi^'^, 
the counterion adsorption in polyelectrolyte solutions was described as an analog of ion- 
pairing in simple electrolyte solutions, and deviations from the Debye-Hiickel limiting laws 
in colligative properties of polyelectrolyte solutions in dilute concentration regime had been 
attributed to counterion adsorption. Ignoring interactions among different chains in di- 
lute solutions, modeling a polyelectrolyte chain by an infinite rod and using bulk dielectric 
constant of the medium for describing electrostatic interactions, counterion adsoprtion (or 
"Manning condensation" ) was predicted to be a result of a singularity in the partition func- 
tion arising due to a singular electrostatic potential near the rod. In particular, it was 
predicted that for monovalent monomers and counterions if Is/b > 1, where Ib is Bjerrum 
length (= e'^/4:7ieoekBT, e being the electronic charge, eo being the dielectric constant of the 
vacuum, e being the dielectric constant of the medium and ksT being the Boltzmann con- 
stant times temperature) and b is the charge spacing on the backbone of the chain, then some 
of the counterions from the solution adsorb on the chain until Is/b = 1. On the other hand, 
if Is/b < 1, the Debye-Hiickel approximation (i.e., the linearization of Poisson-Boltzmann) 
can be used to describe the electrostatics, and there is no adsorption. In other words, coun- 
terion adsorption was described as a kind of phase transition (such as the condensation of 
vapors ) driven by electrostatic energy only. 

Over the years, simulation8^'^'^i^>2i»^i^>^ and experiment3^>^»^>^'^>^»^ have provided an 
insight into the counterion adsorption mechanism in polyelectrolyte solutions. Most impor- 
tantly, it was revealed that a polyelectrolyte chain never attains a perfect rod conformation, 



even in salt-free solutions where the rod conformation is usually expected. Also, it has been 
shown that the effect of ion-pairing^"^ on the backbone and dielectric constantfi^ has to be 
considered in order to treat a realistic polyelectrolyte chain. Furthermore, simulations^*^ 
reveal that the counterions retain their translational degrees of freedom along the backbone 
of the chain after being adsorbed. 

These insights obtained from the simulations and experiments have led to a number of 
theoretical description82Si^i^i22ii2ii2^i2^i21i2^ of the counterion adsorption in polyelectrolyte 
solutions (containing monovalent or multivalent salts) over the last two decades. For an 
isolated, single polyelectrolyte chain, it was shown that there is no singularity in the parti- 
tion function when the polyelectrolyte chain is modeled as a flexible one-- . An important 
prediction of the theory^ emphasized the role played by the dielectric mismatch between the 
local environment of the backbone and the bulk solution in driving counterion adsorption, 
and that had not been highlighted earlier in the literature. In particular, it was shown that 
counterion adsorption arises as an interplay of adsorption energy and translational entropy 
of ions. Entropically, higher degree of ionization, an effect which is opposed by the lower- 
ing of electrostatic energy by the formation of ion pairs at the backbone of the chain, is 
preferred. To carry out the calculations analytically, the Debye-Hiickel potential was used 
to describe the electrostatic interactions between charged species, and a novel variational 
method was used to express the chain conformational entropy and excluded volume effects. 
Finally, the complete free energy of the system (consisting of the chain, the counterions, 
and the solvent) was simultaneously minimized in terms of the degree of counterion ad- 
sorption and the size of the polyelectrolyte chain to self-consistently obtain the equilibrium 
values of the respective variables. The theory is in qualitative agreement with the known 
simulation results. Recently, this single chain theory has been extended to describe the com- 
petitive counterion adsorption^^ phenomenon in the presence of multivalent and monovalent 
counterions. Despite the qualitative agreements to simulations and experiments (especially, 
regarding the phenomena of charge reversal and reentrant transition), the approximations 
used in the theory have not rigorously been assessed so far. 

In this work, we consider the counterion adsorption on a flexible polyelectrolyte chain 
in a spherical cavity in the presence of monovalent salt. Physically, such a situation is 
realized in extremely dilute polyelectrolyte solutions, where inter-chain interactions can be 
ignored safely and a finite volume can be carved out for each chain, and in pores confining 



polyelectrolyte chains, depending on the ratio of the cavity size to the radius of gyration 
of the polymer. We use the self-consistent field theory (SCFT) to compute the equilibrium 
degree of counterion adsorption. SCFT computes the free energy of the system by summing 
over all possible conformations of the chain, and hence, provides a more accurate description 
of the system (in fact, it provides the exact free energy at the mean-field level) compared to 
the variational formalism. Also, as the electrostatics in SCFT is treated at full, non-linear 
Poisson-Boltzmann level, we can assess the validity of the Debye-Hiickel potential to describe 
the electrostatic energy in the variational formalism. Although SCFT provides an accurate 
and clear picture, it is computationally expensive to calculate the degree of ionization due 
to a vast parameter space in the case of polyelectrolytes. On the other hand, the variational 
theory put forward by Muthukumar— is transparent, analytically tractable (to some extent), 
and very inexpensive in terms of the computational needs. Aim of this study is to provide 
a simple, accurate, and easy-to-use method to compute the degree of ionization and assess 
the approximations used in the variational formalism. 

This paper is organized as follows: the theoretical formalisms are presented in Sec. [Til 
calculated results and conclusions are presented in Sec. [Tllland HVj respectively. 

II. COMPARISON OF THEORIES : SCFT AND VARIATIONAL FORMALISM 

We consider a single flexible polyelectrolyte chain of total Kuhn segments, each with 
length /, confined in a spherical cavity of volume Q = AnR^/S. The polyelectrolyte chain 
is represented as a continuous curve of length Nl, and an arc length variable t is used to 
represent any segment along the backbone so that t G [0,A^/]. Now, we assume that the 
chain (negatively charged) is surrounded by ric monovalent counterions (positively charged) 
released by the chain along with ions of species 7 (= — ) coming from added salt so that 
the whole system is globally electroneutral. Let Zj be the valency (with sign) of the charged 
species of type j, and Ug be the number of solvent molecules (satisfying the incompressibility 
constraint after assuming the small ions to be pointlike) present in the cavity. For simplicity, 
we have taken the volume of a solvent molecule (vs) to be equal to the volume of the monomer 
(i.e., Vs = l^). Subscripts p, s, c, + and — are used to represent monomer, solvent , counterion 
from polyelectrolyte, positive and negative salt ions, respectively. 

In order to study counterion adsorption, we use the so-called "two-state" model for the 



counterions so that there are two populations of counterions in the system. One population 
of the counterions is free to enjoy the available volume (called the "free" counterions) and 
the other population is "adsorbed" on the backbone. However, the adsorbed counterions 
are allowed to enjoy translational degrees of freedom along the backbone, maintaining a 
total charge of ef'NZp on the chain, where e is the electronic charge and / is the degree of 
ionization of the chain (i.e., there are —{1 — f)NZp/Zc "adsorbed" counterions on the chain). 
In the literature, this kind of charge distribution has been referred to as a "permuted" charge 
distribution^^. 

For a particular set of parameters, we compute the free energy of the system comprising 
of the single chain, its counterions, the salt ions, and the solvent as a function of the degree 
of ionization (/). As mentioned before, we compute the free energy using two different 
computational frameworks: SOFT— and the variational^^'^^'^^ formalism. In both the 
formalisms, we ignore the electrostatic interactions between solvent molecules and the small 
ions, and model the dielectric constant (e) of the medium to be independent of temperature 
(T) to extract energy and entropy of the system. Also, for comparison purposes, we divide 
the free energy into a mean field part and an additional part, which goes beyond the mean 
field theory. Mean field part is further divided into the contributions coming from the 
"adsorbed" counterions (F*) and the "free" ions, and from the chain entropy etc., (F^)- AH 
contributions are properly identified(or subdivided into) as the enthalpic or(and) entropic 
parts. In all of what follows, the superscript * represents the mean field part. A brief 
description of the derivation for the two formalisms is presented in appendices A and B, and 
the original references^^"^ may be consulted for details. 

To start with, we note that the contributions coming from "adsorbed" counterions are 
the same in both formalisms, and are given by 

F: = E,- TSa, (1) 
Ea = -{l-f)N6lB/l, (2) 
-TSa = iV [/log / + (1 - /) log(l - /)] , (3) 

so that Ea is the electrostatic binding energy of the ion-pairs formed on the polymer back- 
bone due to the adsorption of ions and Sa is the translational entropy of the "adsorbed" 
counterions along the backbone. Parameter 6 = el/eid, reflects the deviation of the dielectric 
constant at the local environment of the chain (e^) from the bulk value (e), and d represents 



the length of the dipole formed due to ion-pairing. 



A. Self-Consistent Field Theory 

Although F* is the same in both formalisms, other contributions involving the "free 
" ions, the chain entropy etc. (i.e., F^) differ from each other significantly in terms of 
computational details. In SCFT, is computed after solving for fields experienced by 
different components in the system, which arise as a result of interactions of a particular 
component with the others. For a single polyelectrolyte chain in a spherical cavity, each 
charged component (monomers and small ions) experiences a dimensionless field (in units 
of ksT/e), which is given by the solution of the Poisson-Boltzmann equation 

Vr^(r) = -47r/BPe(r), (4) 

where Pe(i") = J2j=c + - ^jPji^) + ^vf PiA^) is the local charge density and ppiv) is the 
collective number density of species of type [3 = p, c, +, — . Collective number densities for 
small ions are given by the Boltzmann distribution with the prefactor determined by the 
constraint that the number of small ions are fixed in the system. Explicitly, 



rij exp [—Zjtlj{r)] 



(5) 



J dr exp [—Zjilj{r)] 

for j = c, +, — . In Eq. ( |4]), Pp(r) is the monomer density, which is related to the probability 
of finding a particular monomer (described by the contour variable t) at a particular location 
r, when the starting end of the chain can be anywhere in space (say, q{r, t)), by the relation 
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Furthermore, it can be shown that q{r,t) satisfies the modified diffusion equation^"^ 
dq{r,t) 
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dt 
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g(r,t), tG(0,iV), 
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Wp being the field experienced by monomers due to non-electrostatic interactions. At the 
saddle point, it is given by 



(8) 



where Xps is the dimensionless Flory's chi parameter and ?7(r) is the Lagrange's muhipher 
to enforce the incompressibihty constraint 

+ Ps{r) = Po (9) 

at all points in the system. Here, po is the total number density i.e., po = {N + ns)/Q = l/l^. 
Finally, Ps(r) is the collective number density of solvent molecules given by the Boltzmann 
distribution in terms of the field experienced by a solvent molecule [wg)- Within the saddle- 
point approximation, ps and Wg are related by 

risexp [-Ws{r)] 



Ps[r) 



f dr exp [— Ws(r)] ' 
Xpsl'^pp{r) + ri{r). 



(10) 
(11) 



These equations also close the loop of self-consistent equations, and Eqs. ( 0]- [H]) form the 
set of coupled non-linear equations that defines the system. 

After solving these equations for fields (and in turn, for densities), the free energy at 
the saddle point, F^, is divided into enthalpic contributions due to the excluded volume 
and electrostatic interactions and into entropic contributions because of small ions, solvent 
molecules and the polyelectrolyte chain. Denoting these contributions by E^, Ee, Si, Sg, and 
Sp, respectively, Fj is given by 

F}-Fo = E^ + E,-T{S, + Ss + Sp), (12) 

where Fq = y {Nwpp + ngWss) is the self-energy contribution arising from the excluded 
volume interactions characterized by the excluded volume parameters Wij between species i 
and j. Explicit expressions for different constituents of F^ are presented in Table [T]in terms 
of densities and fields at the saddle point. Within the saddle-point approximation, the total 
free energy (Fscft) of the system is given by Fscft — + F^. In order to compare the 
free energies obtained from SCFT and the variational formalism for a given and R, a 
single Gaussian chain of contour length Nl in the volume fl is chosen as the reference frame, 
whose free energy is taken to be zero. This reference free energy of confinement for a single 
Gaussian chain has been subtracted from the polymer conformational entropy in Table [B 
The free energy of confining a single Gaussian chain with Kuhn segments of length / each 
in a spherical cavity of radius R can be computed exactly and is given by^ 
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TABLE I: Comparison of contributions to in SCFT and variational formalism. 



B. Variational Formalism 

In variational calculations^^, a single polyelectrolyte chain, whose monomers interact 
with the excluded volume and the electrostatic interactions in the presence of the small 
ions is approximated by an effective Gaussian chain, whose conformational statistics are 
dependent on the different kinds of interactions in the system. To compute the equilibrium 
free energy, its variational ansatz is minimized with respect to the variational parameter /i, 
which is related to the radius of gyration [Rg) of the chain by i?^ = Nlli/6. Physically, this 
corresponds to the minimization of the free energy of the single chain system with respect 
to the size of the chain. For the computations of the equilibrium degree of ionization, an 
additional minimization of the free energy with respect to the degree of ionization has to be 
carried out. However, due to the intricate coupling between the size of the chain and the 
degree of ionization, the minimizations have to be carried out self-consistently. 

The variational ansatz of the total free energy (F^ariationai) is given by F^ariationai = 
F* + Fj + AF, where AF involves one-loop fluctuation corrections, addressing the density 
fluctuations of the small ions, to the free energy. As mentioned before, the free energy of the 
adsorbed counterions (i.e., F* ) is the same in both SCFT and the variational formalisms 
(cf. Eq. ([ID). The Fj part of the free energy^ is tabulated in Table [B The function Go (a) 
in Table [T] is a cross-over function given by^"^ 

0o(«) = ^ ( - ^) exp(a)erfc(v/^) + ^ + |- ^-^, (14) 



where a = k'^NIIi/6, and kI is the dimensionless inverse Debye length. Furthermore, h = 
h/h^B = ^b/^- The number of salt ions (n+,n_) are related to the salt concentration (cg) by 
the relation Z^n^ = —Z_n^ = 0.6023csfl, where is in units of moles per liter (molarity). 
Also, all the terms in the free energies are in units of ksT. 

In this work, we have ignored one-loop corrections to the free energy within SCFT. 
However, one-loop corrections to the free energy coming from the density fluctuations of the 
small ions, within the variational formalism, is given by 

iZTT 

where = AnlsifN + n+ + n^)/Q and k is the inverse Debye length. 

C. Numerical Techniques 

We solve SCFT Eqs. [( Hj) - ( [TT]) ] within spherical symmetry (i.e. r —>■ r = |r|), using 
the Dirichlet boundary conditions for q{r,t) and all the fields except rj^r). Also, due to the 
use of spherical symmetry in these calculations, we use 

d^{r) dq{r,t) , ^ f ii + fi(i\ 
— r=o = — ^ r=o= for all t. (16) 

or or 

Starting from an initial guess for fields, new fields and densities are computed after solving 
the modified diffusion and Poisson-Boltzmann equation by finite difference method^. Broy- 
den's method^ has been used to solve the set of non-linear equations. The equilibrium value 
of the degree of ionization (/*) is obtained after minimizing the free energy with respect to 
/. We carry out the numerical minimization of free energy over / using Brent's methocP^. 
The results presented in this paper were obtained by using a grid spacing of Ar = 0.1 and 
contour steps of At = 0.01. 

On the other hand, the self-consistent minimization of the free energy in the variational 
method has been carried out by assuming a uniform expansion of the chain within spherical 
symmetry. In this formalism, the free energy is minimized simultaneously with respect to 
/ and /i, and both these quantities at equilibrium (/*,/*) are computed self-consistently. 
The radius of gyration of the chain (which is confined to a finite volume, Q = AnR^/S) 
is obtained from the equilibrium value of the expansion factor I*. For these calculations, 
the upper bound for the radius of gyration of the chain is specified to be the radius of the 



confining volume (i.e, Rg < R) to mimic the confinement effects. Also, the Kuhn step length 
/ is taken to be unity in both variational as well as SOFT calculations. 

III. RESULTS 

A. Degree of ionization 

We have carried out an exhaustive comparison between the SCFT and variational for- 
malisms by calculating the equilibrium degree of ionization (/*) of a negatively charged 
single flexible polyelectrolyte chain (i.e. Zp = —1, = 1) in the presence of a monovalent 
salt. The equilibrium degree of ionization is determined as a function of the strength of 
the electrostatic interaction (or Coulomb strength) which is proportional to the Bjerrum 
length Ib for a given solvent. For both cases, the effective charge expectedly decreases (Fig. 
[T]) with higher Coulomb strengths that help a progressively larger degree of adsorption of 
counterions on the chain backbone. However, /* obtained from the variational procedure is 
systematically higher than that from SCFT. It is to be noted that the degree of ionization 
is essentially zero in SCFT for experimentally relevant values of Ib/1 (around 3 for aqueous 
solutions), whereas /* is reasonable in the variational theory. Although both theories use 
different approximations and computational procedures, there is one major conceptual input 
that distinguishes these theories. Whereas the variational formalism of Ref. includes the 
density fluctuations of the small ions as one-loop corrections to the free energy, the SCFT 
does not address these fluctuations. In an effort to quantify the consequences of small ion 
density fluctuations and then compare the consequences of the rest of the terms in the 
variational theory against SCFT (which does not contain small ion density fluctuations by 
construction), we subtract AF from Fyariationai and then compute /*. The results are given 
in Figs. [2] and El 

Remarkably, in different conditions corresponding to widely varying degrees of con- 
finement, the /* obtained by the minimization of SCFT free energies is indistinguish- 
able from that obtained using variational free energies without one- loop corrections (i.e., 
^variational — AF). Wc demonstrate this in Fig. [2] where we have plotted /* as a function 
of Ib/1 for different spherical volumes (i.e., different R). Thus we arrive at two conclusions: 
(a) density fluctuations of small ions included in the full variational formalism contribute 



significantly in determining the equilibrium degree of ionization and lead to better values of 
/* than SCFT, and (b) the value of /* is remarkably indistinguishable between the SCFT 
and the variational formalism with deliberate suppression of small ion density fluctuations. 
The first conclusion can be readily rationalized as follows by considering the two curves in 

Fig. m 

The increase in /* with the inclusion of AF can be understood by the fact that the 
density fluctuations of the small ions lower the free energy, and its contribution to the total 
free energy increases with the increase in the number of "free" ions (goes like — n^/^ in salt- 
free case, where n is the number of "free" ions - cf. Eq. ( [T5i) ). For higher values of Ib/^ 
(above 4), all the counterions are adsorbed on the chain so that the degree of ionization of 
the chain is zero irrespective of the density fluctuations. In contrast, for lower values oils /I 
(below 0.5), the chain is fully ionized and the effect of the density fluctuations of the small 
ions on the effective degree of ionization is minimal. However, for the intermediate values 
of the density fluctuations of the small ions affect the degree of ionization significantly 
and non-monotonic deviations from the SCFT results as a function of Ib/1 are observed in 
this regime. Also, term-by-term comparison of the free energy components reveals that the 
discrepancy arises solely due to the term accounting for the density fluctuations of the small 
ions. This disagreement highlights the fact that the effect of density fluctuations of the small 
ions is not included in SCFT within the saddle point approximation. 

The above second conclusion requires further scrutiny. The remarkable agreement be- 
tween the two formalisms is surprising, since these theories use different approximations and 
different computational procedures. In the variational formalism of Ref.^, that is used in 
the present paper, the chain swelling due to electrostatic interaction is assumed to be spher- 
ically uniform at all length scales and at the level of Debye-Hiickel potential between the 
segments. However, this scheme is more tractable analytically with different contributing 
factors (Table [T]) having explicit physical interpretation. On the other hand, in SCFT, the 
electrostatic interaction is at the nonlinear Poisson-Boltzmann level, and the chain expansion 
is addressed at all local length scales through fields generated by intersegment potentials. 
Although the chain conformations are not readily accessible in the standard version of SCFT 
used here, the free energy of the system can be calculated and its resolution into entropic 
and enthalpic parts is possible. In view of such apparently divergent approaches in SCFT 
and the variational formalisms, we now proceed to make quantitative comparisons between 



the two in terms of the various contributing factors. 



B. Term-by-term comparison of free energy: SCFT and variational formalism 

To assess the approximations used in the variational theory and to find out the origin of 
the remarkable agreement in terms of the equilibrium effective charge (/*) obtained from 
SCFT and variational theory (with dehberate suppression of one-loop corrections for small 
ion density fluctuations), we have compared individual contributions to the free energies in 
these two formalisms. Before presenting the numerical results, the role of different contri- 
butions in driving the counterion adsorption can be understood qualitatively by considering 
the following physical picture. 

The driving forces for the counterion adsorption are the formation of ion-pairs due to 
the presence of strong attractive interactions in the process (self-energy of the dipoles) and 
the decrease in intramolecular electrostatic repulsions (compared with a fully ionized chain, 
where these repulsions are the strongest). However, an extensive counterion adsorption on 
the chain backbone is unfavorable due to the loss in translational degrees of freedom of the 
"free" counterions. Another factor, which plays a role in this competition of the energy and 
entropy is the translational entropy of the "adsorbed" counterions. This entropic feature 
alone favors a state where half of the charges on the chain backbone are free and the other half 
adsorbing the counterions. With an increase in the electrostatic interaction strength (i.e., 
Bjerrum length), the driving forces for the counterion adsorption increase and drive more and 
more counterions to the backbone. On the other hand, an extensive counterion adsorption 
leads to the chain contraction due to lower intramolecular electrostatic repulsions (a result 
of the lower number of bare charged sites) even when the electrostatic interaction strength 
is high. However, we show that the counterion adsorption leading to the lowering of the 
effective charge (that decreases the electrostatic energy due to the formation of ion-pairs) has 
a bigger effect than the chain contraction (which affects the chain conformational entropy and 
polymer-solvent entropy) or the increase of intramolecular electrostatic repulsions among the 
unadsorbcd segments as we gradually increase the Bjerrum length. Of course, in addition, 
correlations of small ion density fluctuations also contribute to /*, in the full variational 
calculation. Numerical results on the relative importance of the various contributions to the 
total free energy along with their role in driving the counterion adsorption are presented 



below. 

To start with, in Fig. we have plotted the total free energy calculated in both methods 
for the following set of parameters: Zp = —Z^ = —l,R/l = 10, = 100, Cs = O.IM, 5 = 3 
and Xps = 0.45. It is clear that the total free energies obtained from SCFT and the variational 
theory are in quantitative agreement with each other. 

To analyze this striking agreement between the two methods we focus on the individual 
components of the free energy as tabulated in Table [H In Figs. H] and we have 
compared these different constituents of the free energy obtained from both SCFT and the 
variational formalisms for low monomer densities. It is evident that both theories predict 
that the major contributions to the free energy are from the ion-pair energy (Fig. HH), 
the "adsorbed" counterion translational entropy (Fig. Hb), the polymer- solvent interaction 
energy and the solvent entropy (Fig. Ht); and the "free" ions translational entropy (Fig. 
mi). Contributions due to the chain conformational entropy (Fig. [HK) and the electrostatic 
energy (Fig. [5b) are almost negligible (less than 0.1% in the total free energy) as compared 
to others. For low monomer densities (monomer volume fractions lower than 0.1), the 
dominant contributions to the total free energies come from the polymer- solvent interaction 
energy and the solvent entropy. For the particular single chain dilute system investigated 
here, these contributions account for more than 50% of the total free energy. Although large, 
these contributions are found to be almost insensitive to /. For example, polymer-solvent 
interaction energy changes only by less than O-bhsT when Ib/1 is varied from 0.2 to 5.0. On 
the other hand, /* changes from 1 to in the same range of Ib/1- In fact, the / dependent 
terms which contribute significantly to the total free energy are the ion-pair energy and the 
"free" ions translational entropy. At lower electrostatic interaction strengths (i.e., low Is/l), 
the translational entropy of the "free" ions dominates and at higher electrostatic strengths, 
the ion-pair energy term contributes significantly to the free energy. Together, these two 
contributions account for as high as 99% of the / dependent part in the total free energy (cf. 
Figs. [31 [Hi and [Hi). Relatively very small contributions (~ 1%) to the free energies come 
from the translational entropy of the "adsorbed" ions. We will see below, however, that 
the relative importance of a particular contribution in determining the equilibrium degree 
of ionization is not necessarily related to its actual contribution to the total free energy. 

We now discuss the various trends seen in Figs. [Hand [3 based on conceptual arguments 
aided by the different terms in Table [T] Intuitively, stronger ion-pair energy should promote 



counterion adsorption. As IbII is increased, the energy due to counterion adsorption should 
decrease monotonically, as seen in Fig. H^. From Eq. [21 it is clear that the ion-pair energy 
(note the negative contribution) favors counterion adsorption with a linear dependence on / 
and the Coulomb strength Ib-, hence a progressive gain in adsorption energy with increasing 
Coulomb strength (Fig. H^). On the other hand, the counterion adsorption is opposed by 
the translational entropy of the "free" ions (see the expression for TSi in Table [T]), and hence 
a progressive loss of the part of the free energy related to the translational entropy (Fig. 
mi). The plateau in Fig. HJi arises due to the completion of the adsorption of all counterions 
which limits the loss in the number of free counterions (which is the number of salt ions) with 
increasing Coulomb strength. However, there is no plateau in Fig. H^, because even after 
all counterions are adsorbed, the ion-pair energy continually decreases due to an increasing 
Coulomb strength. In addition, the translational entropy of the "adsorbed" counterions 
{—TSa) drives the adsorption toward /* = 0.5 to optimize this part of the entropy (cf. 
Eq. [3]). Physically, it can be understood from the fact that the complete adsorption of the 
counterions leads to the lowering of the translational entropy of the "adsorbed" counterions 
due to the unavailability of sites. Similarly, a complete desorption of the counterions also 
leads to the lowering of translational entropy of the "adsorbed" counterions due to the 
unavailability of the "adsorbed" counterions on the chain backbone. For a given A^, the 
translational entropy of the "adsorbed" counterions is optimum at /* = 0.5. We note, 
however, that the other two contributions might overwhelm —TSa so that at equilibrium 
it is not necessarily at its minimum (Fig. IHb). With varying Coulomb strength, —TSa 
is minimum at around Ib/1 = 0.8 at which / ^ 0.5, which is prevalently determined by 
the first two components mentioned above. Role of other contributions, i.e., the polymer- 
solvent interaction energies and the solvent entropy (i.e., E^ — TSg), the electrostatic energy 
involving the "free" ions and the monomers {Eg) and the conformational entropy of the chain 
{—TSp) is miniscule in driving the counterion adsorption in a particular direction. However, 
these three contributions dictate the effective size of the chain (through li) at the equilibrium 
(note the dependence of these terms on /i in Table [T]). Further, the equilibrium counterion 
distribution specified by the first three contributions stipulates the actual contributions of 
the last three parts of the free energy at equilibrium. We have noticed before that with 
increasing Coulomb strength (/_b//) the number of free counterions (and, therefore, the 
effective charge of the chain) decreases. Due to a decreasing electrostatic repulsion between 



the monomers, the polymer chain contracts progressively until it reaches its Gaussian size 
at zero effective charge. Consequently, there is less mixing between the polymer and the 
solvent at higher Coulomb strengths leading to a gradual loss of polymer-solvent interaction 
energy (Fig. Hb) which saturates to a plateau when all counterions adsorb, the physical 
condition that creates plateaus in all these curves. Also accompanying the decreasing size 
of the chain there is a gain in conformational entropy (which is maximum at the Gaussian 
size) observed in Fig. EK- In addition, a gradual decrease of the effective charge of the chain 
progressively reduces the electrostatic energy penalty observed in Fig. 03- However, this 
effect is very small compared to the lowering of electrostatic energy due to the formation of 
ion-pairs as mentioned earlier. 

The quantitative agreement between the first three contributions to the free energy in 
two formalisms explains the observed agreement in the results obtained for /* (Fig. [2]). 
Despite having a negligible contribution to the total free energy (less than 0.1%), the 
comparison reveals that the Debye-Hiickel estimate for the electrostatic energy (Ee) used 
in the variational formalism is an overestimation (as large as five times the full, non-linear 
Poisson-Boltzmann at low /_b//). In other words, the Debye-Hiickel approximation underesti- 
mates the degree of screening, which is in agreement with other theoretical^^ and simulation 
results-^. Note that the electrostatic energy in Fig. [Hb includes all the charged species in the 
system except the ion-pairs formed on the chain by the adsorbing counterions. Nevertheless, 
contributions due to the electrostatic energy to total free energy are almost negligible and 
hence, do not affect /* significantly. 

We have also carried out the same comparison between the two formalisms at higher 
monomer densities (above monomer volume fractions of 0.1). It is found that the discrepancy 
in the polymer-solvent interaction energy and the solvent entropy between the two schemes 
is significant (see Fig. All other contributing factors are essentially the same between the 
two theories. The origin of this discrepancy lies in the expansion of the (1 — pp) log(l — pp) 
term, which is carried up to only terms quadratic in polymer density in the variational 
calculations (Appendix-B). The higher order terms in the expansion are ignored in the 
variational calculations to carry out the analysis analytically, which limits the applicability 
of the variational theory to sufficiently low monomer concentrations. The discrepancy clearly 
highlights the breakdown of the variational procedure at high densities and questions the 
use of an effective excluded volume parameter in variational calculations. However, we have 



not made an attempt to compute the boundary of the disagreement between the theories 
because our main focus of this article is on /* which is insensitive to this discrepancy. Also, 
the variational formalism predicts the polymer- solvent interaction energy and the solvent 
entropy to be completely independent of the electrostatic interaction strength Ib/1 in the 
high density regime, in contrast to SCFT predictions of a weak dependence on Ib/1 (see Fig. 
[H]). This is a result of the constraint Rg < R used in variational calculations for mimicking 
the confinement effects and shows the inability of the constraint to capture the confinement 
effects in an appropriate fashion. While the radius of gyration Rg of the chain follows readily 
from li in the variational formalism, it is nontrivial to compute this quantity in SCFT. In 
view of this, we have not addressed Rg in the present paper. 

Finally, we remark on the experimental relevance of the radius parameter R for the 
confining cavity. The variational calculation readily gives Rg without any confinement for 
fixed values of monomer density and other parameters such as Ib, Xps, and Cg. Knowing this 
result, we have investigated the role played by the cavity radius R in the above analysis. If R 
is larger than Rg, then the above conclusions are relevant to unconfined dilute polyelectrolyte 
solutions. On the other hand, if R is less than Rg, then confinement effect is manifest and now 
our results are relevant to a polyelectrolyte chain inside a spherical pore. As an example, for 
N = 100, c, = O.IM, 

Xps — 0.45, and 6 — 3, the calculated value of Rg/l from the variational 
procedure depends on Ib and attains a maximum value of 7.29, whereas R/l = 10.0 in Figs. 
E]- O Therefore, the conclusions drawn above based on these figures are generally valid for 
dilute polyelectrolyte solutions. On the other hand, R/l = 4.0 in Fig. [6], whereas the 
maximum value of Rg would have been 5.92 if confinement were to be absent. Under these 
conditions the conclusions regarding the discrepancy between SCFT and variational theory 
is pertinent to a polyelectrolyte chain confined inside a spherical cavity. 

IV. CONCLUSIONS 

In summary, we have computed the effective charge of a single flexible polyelectrolyte 
chain using SCFT and compared it with the results obtained from a variational theory. It 
is found that for all sets of parameters, the effective degree of ionization (/*) computed 
from SCFT and the variational theory is in quantitative agreement if one-loop fluctuation 
corrections are deliberately suppressed in the latter. The origin of this agreement lies in 



the fact that /* is determined as an interplay of the ion-pair energy and the translational 
entropy of the "adsorbed" counterions as well as of all "free" ions. The conformational 

entropy of the chain, the electrostatic energy involving the "free" ions and the chain, the 
polymer-solvent interaction energy and the solvent entropy do not play significant roles in 
affecting /*. 

The comparison of different components in free energy reveals that the Debye-Hiickel ap- 
proximation underestimates screening effects as compared to the Poisson-Boltzmann theory. 
Despite the fact that there are small discrepancies in the different contributing factors to 
the total free energy, the effective degree of ionization (/*) comes out to be the same in 
SCFT and the truncated variational theory. Furthermore, the density fluctuations of the 
"free" ions, which are included in the full variational theory, are predicted to increase the 
equilibrium degree of ionization. As this latter effect is not captured by SCFT calculations 
within the saddle-point approximation and due to the close agreement between SCFT and 
the variational theory for all other contributing factors, the variational theory appears to be 
a very useful tool for a quick, easy and transparent estimation of /*. 
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Here, we present a brief description of SCFT for a single fiexible polyelectrolyte chain 
having a fixed degree of ionization (= /) in the presence of salt ions and solvent molecules. 
We start from the Edward's Hamiltonian, written as 
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APPENDIX A : SCFT for "permuted" charge distribution 




J]5(pp(r)+p,(r)-po)5Q^ 



Nl 

dty{t)-fNe 



(A-1) 



where R(t) represents the position vector for the t^^ segment and subscripts j = s, c, +, — . 
In Eq. ( I A- II) ksT is the Boltzmann constant times absolute temperature. In writing the 
interaction energies between the polyelectrolyte segments and small ions, we have taken the 
small ions to be point charges so that they have zero excluded volume, and hence, interactions 
are purely electrostatic in nature. As we consider the "permuted" charge distribution, the 
partition function has an additional sum over all possible locations of the "adsorbed" ions 
on the backbone, which appears as an average over the parameter y in Eq. ( lA-ip . We define 
the average over y &s [■ ■ -jy = J dy[- ■ -jgiy), where g{y) = f5{y{t) - 1) + (1 - f)5{y{t)). 
In the above equation, the microscopic densities are defined as 
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where pp{r), pj{r) and Pe(r) stand for the monomers, small molecules (both ions and solvent 
molecules) and the local charge density, respectively. The Dirac delta functions involving 
microscopic densities enforce the incompressibility condition at all points in the system (po 
being the total number density of the system). The delta function involving y is a constraint 
that for all the charge distributions to be considered for one particular value of /, the net 
charge on the chain must be a constant (= fNe). Taking different charge distributions 
of the chain for the same net charge (= fNe) and a particular chain conformation to be 
degenerate, the partition function is divided by the number of ways (p) the "adsorbed" 
counterions can be distributed along the chain. If M out of total sites on the backbone 
are occupied at any particular instance, then p is given by p = A^!/(M!(A^ — M)!) so that 
1 - / = M/N. 

The dimensionless Flory parameter for chemical mismatch, Xps, is given by Xps^^ = 
Wps — {wpp+Wss)/2, where Wpp, Wss, and Wps are the excluded volume parameters, which char- 
acterize the short range excluded volume interactions of type monomer-monomer, solvent- 
solvent and monomer-solvent, respectively. Fq and Ea are the self-energy and ion-pair energy 



contributions, respectively, given by 



— — = Nwpp + UsWss, (A-5) 



Ea 



;i - f)N5lB/l. (A-6) 



where 5 = el/eid, e/ and d being the local dielectric constant and the dipole length, respec- 
tively, is used to characterize the formation of an ion-pair on the bakbone due to "adsorbed" 
counterion. 

Now, using the methods of collective variables and the Hubbard-Stratonovich 
transformational^ for the electrostatic part in Eq. ( lA-ip . the partition function can be 
written as integrals over the collective densities and corresponding fields so that Eq. ( lA-ip 
becomes 



exp (-^^^^) = / D[wp]D[pp]D[i;]D[ri]duD[w,]D[p,] 



Here, Wp, Wg are the fields experienced by the monomers and solvent, respectively, and pp, ps 
represent their respective collective densities. All charged species (excluding the ion-pairs 
formed due to adsorption of counterions) experience a field ip (which is equivalent to the 
electrostatic potential). 77 and u are Lagrange's multipliers corresponding to, respectively, 
the incompressibility and net charge constraints in the partition function. We must stress 
here that this procedure is equivalent to introducing collective fields and densities for small 
ions instead of using the Hubbard-Stratonovich transformation for the electrostatics part. 

Within SCFT, the functional integrals over the fields and the densities are approximated 
by the value of the integrand at the extremum (also known as the saddle-point approxima- 
tion). Extremizing the integrand leads to a number of non-linear equations for the fields and 
the densities. The saddle point approximation with respect to u gives equations similar to a 
"smeared" charge distribution, where every monomer has a charge (= fe). The extremiza- 
tion with respect to ip, Wp, pp, r], Wg, and ps, leads to the saddle point specified by Eqs. ( IH- 
El [HI- [H])- Using these saddle point equations and employing the Stirling's approximation 
for Inn!, we obtain the approximated free energy, i.e., F — Fq ~ ^scf^ presented in section 
III Al after taking ksT = 1. The superscript -k represents the saddle point estimate of the 
free energy. 



APPENDIX B : Variational theory 



In this appendix, we present the procedure to obtain the variational free energy as pre- 
sented in Ref.'^ in the absence of ion-pair correlations. In KefM, it has been assumed that 
the counterions from the polyelectrolyte are indistinguishable from the counterions from the 
salt. So, we start from a partition function similar to Eq. ( IA-1|) with the solvent, counteri- 
ons (from the polyelectrolyte and the salt), colons and the chain as distinguishable species. 
After using collective variables, the partition function can be written as 

exp(-^^^) = y"/^K]/^[pJ/^Wn/^[p,]/^Nexp|-^|, (B-1) 

where j = s, c, — and, where the integral over u has already been evaluated by the saddle 
point method so that the functional h corresponds to a single chain with a "smeared" charge 
distribution. Here, we have introduced collective fields and densities for small ions instead 
of using the Hubbard-Stratonovich transformation for the electrostatic part (as already 
pointed out in Appendix-A). This is the analog of Eq. ( I A- 71) in SCFT. Now, evaluating 
the path integrals over wj by the saddle point method, Eq. ( IB-ll) can be written in 
terms of the densities pj. Functional integrals over t] and ps can be carried out in a trivial 
way. To carry out functional integrals over small ion densities, i.e, the pj log pj terms, which 
emerge after integrations over fields Wj, are expanded up to the quadratic terms after writing 
Pj{r) = nj/Q + 6pj{r) so that / dr6pj{r) = 0, and the resulting integrals are Gaussian. This 
procedure also gives one-loop corrections to the free energy coming from the small ions 
density fluctuations (AF/kBT). Now, expanding the (1 — pp) log(l — pp) term up to the 
quadratic terms in pp, the problem of carrying out the functional integrals over Wp and 
Pp is equivalent to a single chain problem whose monomers interact with each other via a 
renormalized excluded volume parameter and an electrostatic potential. The renormalized 
excluded volume parameter comes out to be = 1 — 2xps and the electrostatic potential 
comes out to be the Debye-Hiickel potential, where the inverse Debye length (k) depends 
on the "free" ions only. Eventually, Eq. ( IB-ll) becomes 

where —TSi is the translational entropy of the "free" ions as presented in Table [T] for the 
variational theory. Now, writing the Hamiltonian of a single polyelectrolyte chain using an 



effective excluded parameter (w) and the Debye-Hiickel potential^^, the functional integrals 
over Wp and pp can be computed using the variational technique developed by Muthukumar— . 
Taking ksT = 1, the results are presented in section IIIBI 
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FIGURE CAPTIONS 



Fig. 1.: Comparison of SCFT and tlie variational formalism (with one- loop corrections) to 
illustrate the effect of correlations among small ions on the effective degree of ionization 
(/*). = -Ze = -1, R/l = 10, iV = 100, Cs = O.IM, Xps = 0.45 and 5 = 3. 

Fig. 2.: Comparison of /* computed using SCFT and the variational formalism (without 
one-loop correction) for different values of R and Ib/1- Zp = —Z^ = —1, N = 100, Cg = 
O.IM, Xps = 0.45 and 6 = 3. Plot for SCFT when R/l = 10 is the same as in Fig. [H 

Fig. 3.: Comparison of total free energies (at equilibrium, for / = /*) obtained from SCFT 
and the variational calculations (without one-loop correction). Zp = —Zc = —1, R/l = 
10, = 100, Cs = O.IM, Xps = 0.45 and 5 = 3. 

Fig. 4.: Comparison of major contributions to the free energies (presented in Fig. [3]) 
obtained from SCFT and the variational formalism, (a) Ion-pair energy contributions 
(Ea), (b) translational entropy of the "adsorbed" counterions {—TSa), (c) polymer- 
solvent interaction energy and solvent entTopY{E^—T Sg); and (d) translational entropy 
of the "free" ions {—TSi). The variational theory captures SCFT results quantitatively. 

Fig. 5.: Comparison of minor contributions to the free energies (presented in Fig. [3]). (a) 
Conformational entropy of the chain {—TSp) and (b) electrostatic energy (Ee). 

Fig. 6.: A discrepancy between SCFT and variational theory for the polymer-solvent in- 
teraction energy and solvent entropy contributions arise at high monomer densities. 
Zp = -Z, = -l,R/l = 4:,N = 100, Cs = O.IM, Xps = 0.45 and 5 = 3. 
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